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ABSTRACT 



Context. During the comparison of synthetic spectra calculated with two NLTE model atmosphere codes, namely TMAP and TLUSTY, 
we encounter systematic differences in the EUV fluxes due to the treatment of level dissolution by pressure ionization. 
Aims. In the case of Sirius B, we demonstrate an uncertainty in modeling the EUV flux reliably in order to challenge theoreticians to 
improve the theory of level dissolution. 

Methods. We calculated synthetic spectra for hot, compact stars using state-of-the-art NLTE model-atmosphere techniques. 
Results. Systematic differences may occur due to a code-specific cutoff frequency of the H I Lyman bound-free opacity. This is the 
case for TMAP and TLUSTY. Both codes predict the same flux level at wavelengths lower than about 1500 A for stars with effective 
temperatures (r cff ) below about 30 000 K only, if the same cutoff frequency is chosen. 

Conclusions. The theory of level dissolution in high-density plasmas, which is available for hydrogen only should be generalized to 
all species. Especially, the cutoff frequencies for the bound-free opacities should be defined in order to make predictions of UV fluxes 
more reliable. 
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1. Introduction 



Stars: atmospheres - Stars: individual: HZ 43, Sirius B - Stars: white dwarfs - Ultraviolet: stars - X-rays: 



NLTE model atmosphere codes for hot, compact stars have 
arrived at a high level of sophistication and are success- 
fully employed for spectral analyses, e.g. the Tubingen 
NLTE Model- Atmosph ere Package TMAlfl dWerner et all2003t 
iRauch & Deetienl | 2003l) in the case of LS V +46°21, the central 
star of Sh 2-216 (kaurii et alj|2007l) . 

In the case of high-gravity stars like LSV+46°21 
(logff = 6.9 [cm/sec 2 !, r e ff = 95 000K) or neutron stars (e.g. 
Suleiman ov & Werner! 120071) . the consideration of the dissolu- 
tion of atomic levels by p lasma perturbation (for details see, e.g., 
lHummer~& Mihalas 1988) is crucial fo r a reliable model atmo- 
sphere calculation dHubenv eT al. 1994). 

iBeuermann et alj d2006l 12008b established the DA-type 
white dwarfs HZ 43 (logg = 7.9, r eff = 51 111 K) and Sirius B 
(log g — 8.6, T e g = 24 897 K) as soft X-ray standards. For a cross- 
calibration between the Chandra LETG + HRC-S, the EUVE 
spectrometer and the ROSAT PSPC, pure hydrogen TMAP 
model atmospheres and synthetic spectra were used. In the case 
of Sirius B with a higher surface gravity g (\ogg — 8.6), the level 
dissolution due to pressure ionization is even more efficient 
(Fig. Q} and thus, it is more important to consider it properly. 



2. Level dissolution 

Recently, Jelle Kaastra (priv. comm.) has drawn our attention 
to a deviation between TMAP and fluxes for Sirius B 
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1 http://astro.uni-tuebingen.de/~rauch/TMAP/TMAP.html 

2 http://nova.astro.umd.edu 
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Fig. 1. Occupation probabilities of the lowest 14 H i levels con- 
sidered in our TMAP model atmosphere calculation for Sirius B. 
Note that the line formation takes place at logra » -5.5 3.0. 



(Fig.|2]l at wavelengths lower than about 1 500A, while in the 
case of HZ 43, the model fluxes are in agreement. 

In the last two decades, we have thoroughly compared TMAP 
and TLUSTY from time to time and found only negligible dif- 
ferences due to different numerical approaches and slightly dif- 
ferent constants. In an investigation of the flux discrepancy, we 
are now able to identify its reason. Both codes, TMAP as well 
as TLUSTY, f ollow a generalized form to consider the level 
dissolution by iHubenv et al.l dl994l) . A hitherto unsolved prob- 
lem, however, is a precise, generalized formulation of the ex- 
trapolation of the hydrogen Lyman bound-free opacity into a 
pseudo-continuum below the unperturbed position of the absorp- 
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Fig. 2. Discrepancy between synthetic spectra for Sirius B calcu- 
lated by TMAP (thick, long cutoff, see text) and TLUSTY (thin, 
short cutoff, Lanz priv. comm.) with the same parameters. 



tion threshold (cf. Dap pen et al.l fl987). TMAP uses a heuristic 
approach with 



yf = f- 
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th,0 



(1) 



where v* and vj h '° are the extrapolated and unperturbed thresh- 
old frequencies, respectively, of level i and n\ is its principal 
quantum number, i.e., n\ — 1 for the Lyman continuum. Since 
this approach results in artificial absorption edges at / = 1 (cor- 
responding to A* = 1215.67. A, "short" cutoff), TMAP intro- 
duced / = 0.5 (corresponding to /t' h = 2431.34 A, "long" cut- 
off) in order to achieve a smooth transition into the continuum 
(Fig.0. 
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Fig. 3. Synthetic fluxes for Sirius B calculated from continuum 
models (no line transitions considered for test reasons) with a 
short (thin line) and a long cutoff (thick). The dashed line shows 
a synthetic spectrum with neglection of the level dissolution. 
Note that the artificial absorption edge at the short cutoff is hid- 
den in the strong L a absorption when those lines are considered 
(dotted). 



TLUSTY uses a different treatment of the continuum. Based 
on the physical picture used to der ive the pseudo-continuum 
formulation (cf. Dappe net al.lll987h . which is valid only near 
the ionization limit, an artificial "short" cutoff of the bound- 
free opacity of the hydrogen Lyman-continuum at = 925 A 
(Lanz priv. comm.) may be chosen (TMAP: = 2431.34 A). 
A test calculation has shown that the model atmosphere fluxes 
of TMAP and TLUSTY agree within 5 % if a long cutoff is used 
by TLUSTY, too (Lanz priv. comm., Fig.|4j. 
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Fig. 4. Comparison of TMAP and TLUSTY (dotted: long cutoff, 
dashed: short cutoff) fluxes in the vicinity of the largest deviation 
(Fig.©. 



Test calculations of pure hydrogen models with TMAP 
via the WWW interface TMAM^] at log g = 8.6 and T eft - be- 
tween 10000 and 50 000 K have shown that deviations be- 
tween short-cutoff and long-cutoff model fluxes are negligible 
at T e ff ~ 30 000 K due to the increasing degree of ionization. This 
explains the good agreement of TMAP and TLUSTY fluxes in 
case of HZ 43. 



3. Conclusions 

Since no reliable theory is available, the choice of different cut- 
off frequencies of the H i Lyman bound-free opacity in the NLTE 
model-atmosphere codes TMAP and TLUSTY demonstrates that 
the estimate of the pseudo-continuum at longer wavelengths is 
presently an uncertainty and definitely deserves further investi- 
gation. However, the necessity of a cutoff in order to avoid an 
unrealistic opacity in the infrared is shown in Fig. [5] 

Moreover, a reliable theory for level dissolution by pressure 
ionization is presently available only for H i. Since this is impor- 
tant for all other species as well, a generalized theory is highly 
desirable. However, this is out of the scope of this work. 
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Fig. 5. Ratios of synthetic fluxes calculated by TMAP. Shown 
are short/long cutoff (dashed), short/extralong cutoff (thin, the 
extralong cutoff frequency is 10 12 Hz), and long/extralong cutoff 
(thick). Note that in the case of the extralong cutoff, the flux ratio 
slightly increases above 100 000 A due to an artificial bound-free 
opacity in the infrared. 



References 

Beuermann, K., Burwitz, V., & Rauch, T. 2006, A&A, 458, 541 

Beuermann, K., Burwitz, V., & Rauch, T. 2008, A&A, submitted 

Dappen, W., Anderson, L., & Mihalas, D. 1987, ApJ, 319, 195 

Hubeny, I., & Lanz, T. 2003, in: Stellar Atmosphere Modeling, eds. I. Hubeny, 

D. Mihalas, K. Werner, The ASP Conference Series Vol. 288, p. 5 1 
Hubeny, I., Hummer, D.G., & Lanz, T. 1994, A&A, 282, 151 
Hummer, D. G., & Mihalas, D. 1988, ApJ, 331, 794 

Rauch, T., & Deetjen, J. L. 2003, in: Stellar Atmosphere Modeling, eds. I. 
Hubeny, D. Mihalas, K. Werner, The ASP Conference Series Vol. 288, p. 
103 

Rauch, T, Ziegler, M., Werner, K., et al. 2007, A&A, 470, 317 
Suleimanov, V, & Werner, K. 2007, A&A, 466, 661 

Werner, K., Dreizler, S., Deetjen, J. L., et al. 2003, in: Stellar Atmosphere 
Modeling, eds. I. Hubeny, D. Mihalas, K. Werner, The ASP Conference Series 
Vol. 288, p. 31 



